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Abstract. The problems and paradoxes of the Lifshitz theory in application to real dielectric 
and semiconductor materials are reviewed. It is shown that the inclusion of drift current of 
conduction electrons into the model of dielectric response results in contradictions with both 
thermodynamics and experimental data of different experimental groups. Physical reasons why 
the problems and paradoxes arise are analyzed and found to be connected with the violation of 
basic applicability condition of the Lifshitz theory, the thermal equilibrium. A recent alternative 
approach to the resolution of the problems based on the inclusion of screening effects and 
diffusion current is considered and demonstrated to be thermodynamically and experimentally 
inconsistent. It is argued that the inclusion of the diffusion current leads to an even deeper 
violation of thermal equilibrium. Phenomenologically, the Lifshitz theory with role of drift 
and diffusion currents neglected is shown to be free of problems and in agreement with both 
thermodynamics and all available experimental data. 



1. Introduction 

It is common knowledge that tlie Lifshitz theory gives a description of the van der Waals 
and Casimir forces between two plane parallel plates (semispaces) in terms of the frequency- 
dependent dielectric permittivities of the plate materials £{uj) [1-3]. Starting in 2000, problems 
and paradoxes with the application of the Lifshitz theory to real materials (at first to metals 
[4] and then to dielectrics [5]) were recognized. Specifically, for metals described by the Drude 
dielectric function a large thermal effect was predicted at short separations [4] in qualitative 
disagreement with the case of ideal metals. This prediction was excluded experimentally at a 
99.9% confidence level [6,7]. At large separations, for the Drude metals one half of the Casimir 
force acting between ideal metals was obtained [4]. The same result was obtained in [8,9] 
by modelling a metal with a set of nonrelativistic hard-core or point-like particles confined by 
walls. In [10] it was shown that the substitution of the Drude dielectric function into the Lifshitz 
formulas results in a violation of the third law of thermodynamics for metals with perfect crystal 
lattices. In [11] it was demonstrated that for metals with impurities the Nernst heat theorem 
is satisfied. This, however, does not solve the problem because in accordance with quantum 
statistical physics the Nernst theorem must be satisfied for metals with perfect crystal lattices 
as well. The situation is even more dramatic for dielectrics, the materials for which the Lifshitz 
theory was originally formulated. It was shown that the inclusion of the dc conductivity arising 
in all dielectrics at nonzero temperature radically changes the theoretical predictions and leads 
to disagreement between the Lifshitz theory and thermodynamics [5,12,13]. 
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In this paper we review recent paradoxial applications of the Lifshitz theory to real materials 
(mostly to dielectrics) and analyze the physical reasons of why the paradoxes appear. We show 
that difficulties arise each time when the drift or diffusion currents are included in the model of 
the dielectric response. In this connection we discuss recent experiments with semiconductors 
[14,15] and dielectric [16] test bodies and show that the measurement results are in disagreement 
with theory taking the dc conductivity into account. Recent theoretical approaches attempting 
to solve problems of the Lifshitz theory by accounting for screening effects and diffusion currents 
[17,18] are demonstrated to be in contradiction with both thermodynamic and experimental 
data. 

Our analysis leads to the conclusion that the roots of the problems are in the violation of the 
applicability conditions of the Lifshitz theory when the role of free charge carriers is included 
in the model of the dielectric response. We argue that in this case one violates the condition of 
thermal equilibrium by admitting a unidirectional flux of heat from the Casimir plates to the 

heat reservoir. 

In Sec. 2 we discuss the consistency of the Lifshitz theory with thermodynamics. Sec. 3 
is devoted to the experimental tests of the Lifshitz theory under different assumptions about 
the model of the dielectric response. Sec. 4 contains the formulation of a phenomenological 
approach of how to apply theory in agreement with both thermodynamics and the experimental 
data. In Sec. 5 we provide thermodynamic and experimental tests for the proposed modification 
of the reflection coefficients due to screening effects and diffusion current. Our discussion and 
conclusions are presented in Sec. 6. 



2. Consistency of the Lifshitz theory with thermodynamics 

The main result of the Lifshitz theory is the expression for the free energy of dispersion 
interaction per unit area of two semispaces at a separation distance a, at temperature T in 
thermal equilibrium with an environment [1-3] 
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Here, k± is the magnitude of a wave vector in the plane of boundary planes (perpendicular to 
the z-axis), prime adds a multiple 1/2 to the term I = 0, = 27rkBTl/h {I = 0, 1, 2, . . .) are 
the Matsubara frequencies, kB is the Boltzmann constant, and the Fresnel reflection coefficients 
along the imaginary frequency axis are given by (n = 1, 2) 
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The free energy of atom-wall interaction, where is the dielectric permittivity of wall 
material, can be expressed as 
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where a(i^;) is the dynamic atomic polarizabihty. 

For dielectrics, the dielectric permittivity along the imaginary frequency axis is commonly 
presented in the form [19] 

where ojj ^ are the oscillator frequencies, gj are the oscillator strength and jj are the relaxation 
parameters. This representation assumes some finite value of the static dielectric permittivity 



eo = e(0) = 1 + V 4 < oo- (6) 



Using this model one can find the asymptotic behavior of the free energy of the Casimir 
interaction between two similar dielectric plates with e^^^ = e^^^ = e at low temperature 
T < Teff = hc/{2akB) [5] 
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T{a,T) « E{a) - ^C(.3)ri{eo + 1) (^^J . (7) 

Here, E{a) is the Casimir energy at zero temperature, C{z) is the Riemann zeta function, and 
ro = (eo ~ l)/(£o + 1) is the value of the TM reflection coefficient at zero frequency. A similar 
asymptotic expression can be obtained from Q for the free energy of the atom-wall interaction 
[20] 

T^(a,T)^E^{a)-^^^aiO)CD{eo)[^jr^) , (8) 

where Cd{£:o) is some constant found in [20]. 

Equations ([7|) and ([8]) allow one to find the asymptotic behavior of the Casimir entropy in 
the plate-plate and atom-plate configurations 

at low temperature. Thus, from ([7]) and ^ one obtains 

As is seen from (llOh . in both cases the Casimir entropy goes to zero when T vanishes, i.e., the 
Lifshitz theory with the dielectric permittivity dS]) satisfies the Nernst heat theorem. 

The model of the dielectric response ([5]) is, however, not exact. It does not take into account 
that at T 7^ all dielectric materials possess small but physically real dc conductivity cr{T). For 
all dielectric materials, the general behavior of ct(T) with vanishing temperature is given by 

cr(r) ~ exp f--^^ when T ^ 0, (11) 
V kbTJ 

where the physical meaning of the coefficient C is different for different classes of dielectrics 
[21,22]. The dc conductivity of dielectrics is usually taken into account by means of the Drude- 
like term (see, e.g., [23]) 

m,T)=em + '^. (12) 



It would be reasonable to expect that the substitution of the dielectric permittivity (I12p . 
instead of ([5]), into the Lifshitz formulas ([T]) and (jl]) leads to only minor corrections to the 
asymptotic results (j7]) and ([8]). In reality, however, this substitution radically changes the 
asymptotic behavior of the free energy at low temperature: 
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Here, Lin{z) is the polylogarithm function, T{a,T) and !F^{a,T) are given by ([7]) and 
respectively, and RiT), R{T) are exponentially vanishing functions when T ^ 0. In this case the 
definition of the Casimir entropy ([9|) results in the following nonzero values at zero temperature 
[5,20] 
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This is in violation of the third law of thermodynamics (the Nernst theorem) because not only are 
these values not equal to zero, but they also depend on the parameters of the system (separation 
distance, static dielectric permittivity and atomic polarizability) . 

The obtained result can be considered as a puzzle because instead of more accurate asymptotic 
expressions one arrives to the conclusion of an apparent contradiction of the Lifshitz theory with 
the fundamental principle of thermodynamics. 



3. Experimental tests 

Rapid progress in the measurements of the Casimir force is very helpful for the resolution of 
the theoretical controversies. As mentioned in the Introduction, the experiment [6,7] excluded 
the use of the Drude model to describe the thermal Casimir force between metal plates. It is 
important to note, however, that a much smaller thermal effect, as predicted by the generalized 
plasma model [24,25], was not registered in that experiment. Here we consider the results of 
two other experiments involving semiconductor and dielectric materials [14-16]. Both of them 
are accurate enough to resolve between the theoretical approaches including and neglecting 
the dielectric dc conductivity. What is more, one of them [16] provides the first experimental 
observation of the thermal effects in Casimir physics. 

We start with the observation of modulation of the Casimir force between a gold coated 
sphere and a Si plate with laser light [14,15]. In this experiment, a B doped p-type Si plate 
was illuminated with laser pulses and the difference of the Casimir forces between a sphere 
and a plate in the presence and in the absence of light has been measured. In the absence 
of laser pulse the concentration of charge carriers in the Si plate was n = 5 x 10^^ cm"^ and 
in the presence of light ni = 2.1 x 10^^ cm"^ or 71-2 = 1.4 x 10^^ cm"^ for the two different 
absorbed powers P^i = 9.3 mW and = 4.7 mW. The directly measured quantity was, thus, 
AF{a) = F^{a) — F{a), where and F are the Casimir forces in the presence and in the 
absence of laser light on the plate, respectively. The experimental data were compared with 
the Lifshitz theory either neglecting or including dc conductivity of the high-resistivity Si in the 
dark phase. 

In Fig. l(a,b) the mean experimental data {AF(a)) are shown as dots for the absorbed powers 
Pmi and Pw2, respectively. The calculational results for AF{a), using the Lifshitz theory with 
neglect of the dc conductivity of the high-resistivity Si (dark phase), are shown as the solid 
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Figure 1. Differences of the Casimir forces in the presence and in the absence of hght versus 
separation for different absorbed powers: (a) 9.3 mW and (b) 4.7 mW. The mean measured 
differences are shown as dots. Differences calculated using the Lifshitz formula at T = 300 K 
and T = are shown as the solid and short-dashed lines, respectively, and those calculated 
including the dc conductivity of high-resistivity Si as the long-dashed lines. 



lines. Similar calculational results using the Lifshitz theory with included dc conductivity ()12p 
in the dark phase are shown as the long-dashed lines. As is seen in Fig. l(a,b) the solid lines are 
consistent with the data whereas the long-dashed lines are excluded over the wide separation 
region. In the same figure, the short-dashed lines show the computational results obtained from 
the Lifshitz formula at zero temperature. As is seen in the figure, the computational results at 
T = are consistent with the data. These results do not depend on whether the dc conductivity 
was included or neglected. The experiment under consideration is precise enough to exclude the 
large thermal effect, as obtained when the dc conductivity in the dark phase is included, but is 
not of sufficient precision to measure the small thermal effect obtained in the case when the dc 
conductivity in the dark phase is neglected. 

For illustrative purpose, in Fig. 2 we demonstrate the comparison between experiment and 
theory over a more narrow separation interval from 100 to 150 nm where only each third 
experimental point is plotted together with its error bars [itAa, ibA*°*(AF'^^P*)] shown as crosses 
(the absorbed power is 8.5 mW). Note that these errors are entirely experimental. They are 
determined at a 95% confidence level as a combination of random and systematic errors and are 
not connected with any theory. The theoretical force differences AF^^'^"^ computed at T = 300 K 
with neglected and included dc conductivity in the dark phase are shown as the solid and dashed 
lines, respectively. From Fig. 2 one concludes that inclusion of the dc conductivity in the dark 
phase is excluded by the experimental data at a 95% confidence level within the separation 
region from 100 to 150 nm. 

Note that computations of the Casimir force between an Au coated sphere and Si plate were 
done using the proximity force approximation [26] 

F{a,T) =2TrRJ'{a,T), (15) 

where !F{a,T) is given by ([T|). The dielectric permittivities ^'^"'^^(i^) and e''^)(i^) for Si and Au, 
respectively, were obtained from the tabulated optical data [23] by means of the Kramers-Kronig 
relations. Practically the same results can be obtained using the analytical representations for 
e(")(iC). For Au high precision six-oscillator model is considered in [7]. For Si, the imaginary 
part of the dielectric permittivity can be modeled as a constant within some frequency region 
[27] 

ImsW(^) = (^" ^o<^<a;i. 
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Figure 2. Mean experimental differences of the Casimir forces with their experimental errors 
are shown as crosses. Solid and dashed lines represent the theoretical differences computed 
at T = 300 K using the model of dielectric permittivity of high-resistivity Si neglecting and 
including the dc conductivity, respectively. 

where ujq = 3.22 eV, wi = 4.62 eV, e = 48. Using the Kramers-Kronig relation, this results in 



The use of the approximation (|17p instead of the tabulated optical data [23] leads to less than 
1% deviations in the Casimir force (jlSp between an Au sphere and Si plate. At separations 
a = 50, 100, 200, 400 and 600 nm this deviation is equal to -0.58%, 0.30%, 0.70%, 0.84% and 
0.86%, respectively. 

Another analytic formula representing the dielectric permittivity of Si along the imaginary 
frequency axis takes the form [28] 



where ujq = 6.6 x lO^'^^rad/s, ^'•^^(iO) = 11.87 and e^"'^^(ioo) = 1.035. This representation also 
leads to less than 1% differences in the case of Au sphere above Si plate in comparison to 
computations using the tabulated optical data for Si. Thus, at separations a = 50, 100, 200, 
400 and 600 nm the relative deviations between both computations are equal to 0.75%, 0.74%, 
0.62%, 0.55% and 0.52%. Because of this, both the tabulated optical data for Si and different 
analytical representations for its dielectric permittivity can be used in the computations within 
1% precision. 

Another important experimental test is the measurement of the thermal Casimir-Polder 
force through center-of-mass oscillations of the trapped Bose-Einstein condensate [16]. In this 
experiment, not only were different models of the dielectric response compared with the data 
[16,29], but the thermal effect was measured for the first time in Casimir physics. The condensate 
of Rb atoms was at a distance from 7 to 11 /xm from a fused Si02 dielectric substrate. In the 
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original publication this substrate was considered as a dielectric material with the static dielectric 
permittivity £^^^(0) = 3.81. However, at room temperature (in an equilibrium situation the 
temperature of the substrate and environment were equal to Tg = Tg = 310 K) Si02 has 
nonzero conductivity which is ionic in nature and varies from 10~^ s~^ to 10^ s~^ depending 
on the concentration of alkali ions which are always present as trace constituents [30,31]. This 
conductivity was taken into account in [29] in accordance with ()12p . 

The directly measured quantity in the experiment [16] was the relative shift of the frequency 
of center-of-mass oscillations of the Bose-Einstein condensate due to the presence of the Casimir- 
Polder force 

ujQ 2mujQ 

Here, m = 1.443 x 10~^^ kg is the mass of the Rb atom, loq is the unperturbed trap frequency 
ujQ = 2-7T X 229rad/s, F^{z,T) in the equilibrium situation is given by the negative derivative of 
J-^i^z, T) defined in ^ with respect to z, and the averaging is done over the period of condensate 
oscillations and over the volume of the condensate cloud. 

The experimental data for the relative frequency shift 7^ are presented as crosses in Fig. 3. 
The arms of the crosses show the total experimental errors determined in [16] at a 67% confidence 
level individually for each data point. The results of theoretical computations using (|19p with 
neglected dc conductivity of Si02 substrate are shown by the solid line [16,29]. As is seen in the 
figure, theory neglecting the dc conductivity of fused silica is in a very good agreement with the 
data. In the same figure, the results of theoretical computations by (llOp with dc conductivity 
of Si02 included are shown by the dashed line [29] . As is seen in the figure, the dashed line is 
in evident disagreement with the first two experimental points. 

The experiment was repeated for two more times in situations out of thermal equilibrium 
{Ts = 479 K, Te = 310 K and Ts = 605 K, Te = 310 K), i.e., with the temperature of the plate 
higher than the environment temperature. In this case the Casimir-Polder force F^{z,T) in 
(fT9]) should be replaced with 

F^{z,Ts,Te) = F^{z,Te) + F-^(z,r5) - F'''^z,Te), (20) 

where the explicit form of the function F^'^'^{z,T) can be found in [16,29,32]. 

In Fig. 4 the experimental data for 7^ are shown for (a) Ts = 479 K and (b) Ts = 605 K as 
crosses with their arms indicating the experimental errors. As in Fig. 3, the theoretical results 
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Figure 3. Fractional change in the trap frequency versus separation in thermal equilibrium 
with Ts = Te = 310 K computed by neglecting (solid line) and including (dashed line) the 
conductivity of the dielectric substrate. The experimental data are shown as crosses. 
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Figure 4. Fractional change in the trap frequency versus separation out of thermal equilibrium 
with (a) Ts = 479 K, Te = 310 K and (b) Ts = 605 K, Te = 310 K. Computations are done by 
neglecting (solid line) and including (dashed line) the conductivity of the dielectric substrate. 
The experimental data are shown as crosses. 



computed with the dc conductivity of Si02 plate neglected and included are shown by the 
solid and dashed lines, respectively. As is seen in Fig. 4(a), the dashed line is in disagreement 
with the three experimental points, and the other two only touch it. In Fig. 4(b) all data 
points exclude the theoretical approach represented by the dashed line. Thus, the inclusion of 
the dc conductivity of the dielectric substrate into the Lifshitz theory leads to disagreement 
with the experimental data. The calculational results using the Lifshitz theory with the dc 
conductivity neglected are in very good agreement with the experimental data. We emphasize 
that the theoretical results obtained with dc conductivity included do not depend on the value 
of conductivity, but only of the fact that it is nonzero. Importantly, the inclusion of conductivity 
does not influence the nonequilibrium contributions to jz- Thus, the puzzling effect comes only 
from the contribution to jz given by the Lifshitz theory. 



4. Phenomenological approach and its justification 

In the above, both theoretical and experimental problems are elucidated arising in the Lifshitz 
theory when applied to dielectric materials at nonzero temperature with dc conductivitincludedy. 
Recently, these problems were discussed also in [33] but no solution was found. In [33] only 
"weakly conducting" materials were considered which possess nonzero conductivity at zero 
temperature. These are in fact metallic-type semiconductors with a dopant concentration above 
critical. To satisfy the Nernst heat theorem, it was suggested that there is nonzero residual 
relaxation at T = 0, i.e., the crystal lattice is not perfect. Materials of such type have nothing 
to do with dielectrics whose conductivity vanishes when temperature goes to zero. The formalism 
of [33] is also in contradiction with precision experiments dealing with metallic test bodies [6,7]. 

To avoid contradictions of the Lifshitz theory with the principles of thermodynamics and 
experimental results, one can propose the following phenomenological prescription. In the 
applications of the Lifshitz theory to all materials having zero conductivity at zero temperature, 
the presence of free charge carriers at nonzero temperature should be disregarded. The dielectric 
permittivity of such materials should be described by equation ([5|). In the applications of the 
Lifshitz theory to all materials having nonzero conductivity at zero temperature, free charge 
carriers should be included in the framework of the generalized plasma-like model 

where all ujj ^ 0. This prescription is consistent with thermodynamics and all available 
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experimental data [22,24,25,34]. 

The physical reasons of why the Lifshitz theory does not allow the inclusion of real 
conductivity processes can be understood as follows. Lifshitz derived his famous formulas 
under the condition of thermal equilibrium. This means that not only T = const, but also 
all irreversible processes connected with the dissipation of energy into heat have already been 
terminated [35,36]. The Drude-like dielectric function (I12p is derived from the Maxwell equations 
with a real drift current of conduction electrons j = gqE initiated by the external electric 
field E [24]. The drift current is an irreversible process which takes a system out of thermal 
equilibrium. This current leads to Joule heating of the Casimir plates (Ohmic losses) [37]. To 
preserve the temperature constant, one should admit that there exists an unidirectional flux of 
heat from the medium to the heat reservoir [38]. Such interactions between a system and a heat 
reservoir are prohibited by the definition of thermal equilibrium [39]. One can conclude that 
the substitution of the Drude-like dielectric permittivities into the Lifshitz formulas violate the 
applicability conditions under which this theory was developed. Then it is not surprising that 
such substitution results in problems with thermodynamics and contradictions with experimental 
data. From this point of view, the Lifshitz theory is free of problems and paradoxes which arise 
only if it is applied incorrectly. 

Another question to be answered is the following. One could agree with the fact that the 
conductivity processes in real materials connected with the drift current of conduction electrons 
violate the applicability conditions of the Lifshitz theory and, thus, are not described by this 
theory. These processes, however, exist in reality. Because of this should they be taken 
into account in some future general theory of dispersion forces? At the moment we have no 
answer to this question. All computations on the basis of the Lifshitz theory for dielectrics 
with dc conductivity neglected and for metals by using the plasma-like dielectric permittivity 
()2ip are in agreement with both thermodynamics and experimental data. This suggests that 
the conductivity properties are not related to dispersion forces. More precise experiments will 
probably bring additional information to bear on this subject. 



5. Approaches based on the modification of reflection coefficients due to screening 
effects 

Recently it was proposed [17,18] that the problems of the Lifshitz theory described above may 
have an alternative solution if the screening effects due to free charge carriers are taken into 
account. This approach includes both the drift and diffusion currents of free charge carriers 
through use of the transport Boltzmann equation [18]. It uses the standard Lifshitz formulas 
([1]) and (j4]) for the free energy of wall-wall and atom-wall interaction with the TE reflection 
coefficient, fxE) as defined in ([2]) within the Drude model approach, but with the modified TM 
reflection coefficient 

rM\a -k-^ em-em 
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where e(i^) is the permittivity of the bound core electrons given in ([5]). The quantity r]{iS,) is 
defined as 
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where 1/k is the screening length. If the charge carriers of density n are described by the classical 
Maxwell-Boltzmann statistics, as assumed in [18], one gets the Debye-Hiickel screening length 



[40,41] 

Assuming Fermi-Dirac statistics, one arrives at the Thomas-Fermi screening length [40,41] 

R^p=^=(£^y^\ (26) 

Here, e is the electron charge, Eq = e(0) is the dielectric constant due to the core electrons, and 
Ep = hujp is the Fermi energy. 

Note that the paper [18] applies this approach to intrinsic semiconductor media. Bearing 
in mind, however, that the classical Drude dielectric permittivity is exploited, it is justified to 
apply it for metals and doped semiconductors (both dielectric and metallic type) as well, by 
substituting the respective screening length defined by ([25]) or (p6]) . 

Now we consider whether or not the approach taking the screening effects into account is 
consistent with thermodynamics. We concentrate our attention on the case of dielectric materials 
whose conductivity goes to zero when the temperature vanishes. We start from the case of two 
similar dielectric semispaces [e^^^ = e^"^^ = e]. The Casimir free energy per unit area is given by 
([T]) with the reflection coefficients r!^^ = r!^^ = tte from ([25]) . and r!^y^ = r!^y^ = ttm from 



The derivation of the asymptotic behavior of the Casimir free energy and entropy at low 
temperature is performed in perfect analogy with the derivation of ([7]) and p3|) in [5]. First we 
introduce the small parameter (^ > .^i) 

. 47ro-(r) , , 
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which goes to zero due to (llip when the temperature vanishes. Then we expand the reflection 
coefficients rTM,TE at all nonzero Matsubara frequencies in powers of the parameter (3 



v/y^ + (e-l)C^[ey + ^y^ + {e-l)e] 



Vr + (e - i)C^[y + vr + (e - i)C^] 

Here, the reflection coefficients rxM.TE are defined by ([2]) with the dielectric permittivity ([5]). 
They lead to the asymptotic behavior of the Casimir free energy per unit area given in ([7]). The 
dimensionless variables C and y are connected with other dimensional ones by the equations 

y = 2aq, ( = - = —, (29) 
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where q = q{S,) is defined in ([3]). 

Now we substitute ([28p in ([1]) and arrive at the following expression for the Casimir free 
energy per unit area taking the screening effects into account 

^(a,r) = Tia,T) + {/"^^^/^^ " ^oiy)^~'] + ^'^(^o) + QiT)] , (30) 



where 
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and Q{T) vanishes exponentially when T — > 0. The free energy T{a,T) in (I30p is given by ([7]), 
and K in (fSTI) is defined in (f25]) . Calculating the negative derivative of both sides of ([30]) with 
respect to T, we obtain the asymptotic behavior of the Casimir entropy at low temperature 
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where S{a,T) is defined in ()10p . It can be easily seen that the last three terms in the curly 
brackets on the right-hand side of this equation go to zero when T goes to zero for any dielectric 
material. 

The behavior of the first two terms in the curly brackets on the right-hand side of (|32p when 
T goes to zero is more involved. If n(T) exponentially decays to zero with vanishing temperature 
(as is true for pure insulators and intrinsic semiconductors), then according to ([^^ so does k. As 
a result, fo(?/) — >• ro = {eq — 1)/ (^o + 1) and the first two terms in the curly brackets cancel. Then 
the Casimir entropy S{a,T) goes to zero when T vanishes following the same law as S{a,T), 
i.e., as T^. This means that for insulators and intrinsic semiconductors the formalism under 
consideration is in agreement with the Nernst heat theorem. 

However, there is a wide class of dielectric materials (such as doped semiconductors with 
dopant concentration below critical and solids with ionic conductivity) for which n does not go 
to zero when T goes to zero. In fact, conductivity can be presented as [41] 

cr{T) = n |e| (33) 

where /i is the mobility of charge carriers. Although o"(T) goes to zero exponentially fast for all 
dielectrics when T goes to zero, for most of them this happens due to the vanishing mobility. 
For instance, as mentioned in Sec. 3, the conductivity of Si02 used as the plate material in the 
experiment [16] is ionic in nature and is determined by the concentration of impuritues. For all 
such materials, in accordance with psp . k — > oo when T — > 0. As a result, ro(y) 1 when T 
goes to zero in accordance with (|3T]) . In this case we obtain from (|32|) 



^K0) = ^[C(3)-Li3(rg)] >0 (34) 

in violation of the Nernst heat theorem. This means that for a wide class of dielectric 
materials the proposed approach taking the screening effects into account is thermodynamically 
inconsistent in the same way as the standard Lifshitz theory with the dc conductivity included 
[compare with (I14p ]. 

We emphasize that the existence of dielectric materials for which n does not go to zero but 
/X does go to zero when T vanishes demonstrates that the reflection coefficient ()22p at ^ = is 
ambiguous. In reality, for such materials rTM(0, A;_l) — > 1 when T and /U simultaneously vanish. 
This is because k — > oo when T — > in disagreement with physical intuition that there should 
be no screening at zero mobility. This ambiguity is connected with the break of continuity of 
the reflection coefficient rTM(iCi^±) &t the point ^ = 0, T = 0. If one takes the limit T ^ 
first, keeping ^ = const ^ 0, the standard Fresnel reflection coefficients ttm from ^ with no 
screening is reproduced. This property is preserved in the subsequent limiting transition ^ — > 0. 
As was noted in [42] on general grounds, the violation of the Nernst heat theorem is caused by 
the breaks in the continuity of the reflection coefficients at the origin of the (^, T)-plane. 



The theoretical approach taking the screening effects into account can also be applied to 
describe the Casimir-Polder atom-wall interaction. In this case it is also inconsistent with 
thermodynamics for dielectric materials with a nonvanishing density of charge carriers. To find 
the asymptotic behavior of the free energy at low temperature, one should substitute into ([!]) 



the reflection coefficients r, 



(1) 

TM,TE 



^TM,TE, as deflned in (f28ll . The result is 



ksTaiO) 
8a3 



y^dyroiy)e-y -2ro + QiiT) 



(35) 



where the asymptotic behavior of J-^ia, T) is presented in ([S]), and Qi{T) vanishes exponentially 
when T — > 0. Calculating the negative derivative of both sides of (|35p with respect to T, one 
obtains 



S^{a,T) = S^{a,T) + 



8a3 



y^ dy ro(y)e" 



2rn 



(36) 



dy 



dT Jo "VFT(2H^(eoVFTT2^ + y)^ 



+ Qi(r) + rQ;(r) 



where S^{a,T) is defined in (jlOp . For dielectric materials with exponentially decaying n(T) in 
the limit of low T, fo{y) —>■ tq and all terms on the right-hand side of ([36]) added to S'^{a,T) 
go to zero when T —)■ 0. As a result, S'^{a,T) goes to zero as and the Nernst theorem is 
satisfied. On the contrary, for all dielectric materials with nonvanishing n we get ro(y) ~^ 1 
when T —f and in this case (1361) results in 



5^(a,0) 



kBa{0) 
4a3 



(1 - ro) > 0. 



(37) 



This means that the Nernst theorem is violated in the same way as in (1140 for the standard 
Lifshitz theory with the dc conductivity included. 

The theoretical predictions following from the approach proposed in [18] can be compared 
with the measurement data of the experiment on the modulation of the Casimir force with light 
(see Sec. 3 for details). In Fig. 5(a) and 5(b) we present the experimental data of this experiment 
obtained for the absorbed powers P^ji and Pw2, respectively, with the experimental errors of force 
measurements determined at 70% confidence level. The theoretical difference Casimir force, 
computed using the standard Lifshitz theory with the dc conductivity of high-resistivity Si in 
the dark phase neglected, is shown by the solid lines. The theoretical difference Casimir force 
computed by including the screening effects at all Matsubara frequencies, as described above, is 
shown by the dashed lines. It can be seen that the standard Lifshitz theory is consistent with 
data, whereas the theoretical approach taking the screening effects into account is excluded by 
the data at a 70% confidence level. 

In Sec. 4 the contradictions with thermodynamics and the experimental data of the standard 
Lifshitz theory with the dc conductivity included are explained based on the violation of its 
applicability condition, i.e., of thermal equilibrium, in the presence of drift current. In a 
similar way, the approach taking the screening effects into account applies the Lifshitz theory to 
physical phenomena involving both the drift and diffusion currents. The latter is caused by the 
nonequilibrium distribution of charge carriers in an external field, i.e., by a physical situation 
out of thermal equilibrium. This is because the diffusion current is determined by a nonzero 
gradient of charge carrier density, whereas for homogeneous systems in thermal equilibrium the 
charge carrier density must be homogeneous. Thus, the problems and paradoxes arising in this 
case are also explained by the use of the Lifshitz theory outside of its defined application region. 
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Figure 5. Differences of the Casimir forces between an Au-coated spliere and a Si plate, in the 
presence and in the absence of laser light on the plate, versus separation for the absorbed power 
of (a) 9.3 mW and (b) 4.7 mW. The experimental data are shown as crosses. Solid and dashed 
lines indicate the theoretical results calculated using the standard Lifshitz theory with the dc 
conductivity of Si in the dark phase neglected, and the generalization of this theory [17,18], 
respectively. 



An earlier attempt to incorporate the screening effects in the theory of dispersion forces was 
undertaken in [17]. In this approach the effect of screening was taken into account only for 
the static field in the interaction of an atom with a dielectric wall. Bearing in mind that in 
accordance with Q the TE reflection coefficient at zero frequency does not contribute to the 
result, the single modified coefficient was [17] 



r^^^iO, k^) = ^1 ^ ^ , (38) 

Box k] + + k^ 



where k is given in ()25p . It is easily seen that (I38p coincides with the general TM reflection 
coefficient fTM(iC) k±) defined in (f22|) if one puts ^ = 0. 

The reflection coefficient ()38p can be obtained by the formal replacement of the dielectric 
permittivities depending only on the frequency with the permittivities depending both on 
frequency and the wave vector in the Fresnel reflection coefficient of an uniaxial crystal. To 
confirm that this is the case, one should introduce nonisotropic permittivities £x = £y £z and 
replace the coefficients ttm and tte in (12]) with 



^/ex{-^iyz{ii)q + kz q + kx 

Here, k^ and kx are defined by Q with the replacement of e(iC) by £zO-£,) and ex(iO) respectively. 
For Ex and Ez depending only on the frequency, ([39]) is commonly used for the description of 
uniaxial crystals (see, e.g., [43,44]). Then ([38]) follows from ([39]) if one puts 



ExiO) = Eo, e.(0) = e,(0, k^) = eq [ I + ] , (40) 




where now ez(0) depends on the wave vector, i.e., spatial dispersion is present. Thus, to 
obtain ()38p one substitutes into the Presnel reflection coefficient (I39p derived from the standard 
continuity boundary conditions a dielectric permittivity depending on the wave vector (see 
[45,46] why this is inappropriate). 

We emphasize that the general reflection coefficients taking the screening effects into account 
[the transverse electric from (j2]) with the dielectric permittivity ()23p and the transverse magnetic 
from (j22p ] cannot be represented in the form of (j39p . To prove this statement, we note that from 
the equality tte = r^E^ it follows that e(i^) = £xi^O- Thus, in accordance with ([^5]) . ex(iO ~^ oo 
when ^ ^ 0. Keeping in mind that £2(0) > 0, from (j39p one obtains rjl^{0,k±) = 1. This 
is in contradiction with the above result that at zero frequency the reflection coefficient (|22p 
coincides with the reflection coefficient (|38p . This means that the general approach taking the 
screening effects into account [18] cannot be described in terms of any dielectric permittivity 
which includes spatial dispersion. 

The approach based on the modification of only one reflection coefficient (|38p meets the same 
thermodynamical difficulties as the general approach considered above. Thus, if one computes 
the Casimir entropy with the reflection coefficient (j38p and all the other standard ones, as 
given by ([2]), it goes to zero when the temperature vanishes for dielectrics with exponentially 
vanishing charge carrier density n. However, for dielectrics with nonvanishing n the result (j37p 
is reproduced, i.e., the Nernst heat theorem is violated [20]. 

The modiflcation of the TM reflection coefficient in accordance with (j38p can be verified 
experimentally by using the measurement data of the optical modulation experiment discussed 
in Sec. 3. In theoretical computations at zero frequency the expression (i38]l was used. At 
all nonzero Matsubara frequencies the standard reflection coefficients ([2]) were substituted into 
the Lifshitz formula. The computational results using the standard reflection coefficients at 
all nonzero Matsubara frequencies coincide with those obtained with the general reflection 
coefficients r^TM(iCj^±) and r^TM(i?5^±)- Then we return to the dashed lines in Fig. 5(a,b) 
representing the difference Casimir force computed with the inclusion of screening effects. More 
generally, it can be easily seen that for materials with sufficiently small charge carrier density 
(for instance for Si with n < 10^^ cm~^) the general reflection coefficient ()22p including the 
screening effects leads to approximately the same values at room temperature, as the standard 
one at all nonzero Matsubara frequencies. Thus, the theoretical results computed by using the 
modified TM reflection coefficient (|38p are excluded by data at a 70% confidence level. 

6. Conclusions and discussion 

In the above we have considered the application of the Lifshitz theory to real metals, dielectrics 
and semiconductors. When the conductivity processes, connected with the presence of free 
charge carriers, are included into the model of the dielectric response, problems and paradoxes 
in the application of the Lifshitz theory are immediately apparent. They are manifested in the 
contradictions of the Lifshitz theory with both thermodynamics and experiment. We have listed 
several measurements performed by different experimental groups with metals, dielectrics and 
semiconductors which are in contradiction with the theoretical results taking into account the 
drift current of conduction electrons. 

In searching for reasons of why the problems and paradoxes arise, it was emphasized that 
the Lifshitz theory was formulated for systems in thermal equilibrium. As a result, this theory 
is not applicable to physical systems where processes go on which violate thermal equilibrium. 
Speciflcally, in thermal equilibrium all processes of heat transfer, diffusion and others connected 
with the dissipation of energy are terminated. A question then arises if the drift current of 
charge carriers which leads to a unidirectional flux of Joule heat from the Casimir plates to the 
heat reservoir is compatible with the state of thermal equilibrium. Our conclusion is that it is 
not compatible with thermal equilibrium and, thus, cannot be considered in the framework of 



the Lifshitz theory. The inclusion of conductivity properties in the Lifshitz theory is in violation 
of its basic applicability condition. 

We have analyzed recent attempts to solve problems and paradoxes using screening effects 
and diffusion processes. We have shown that these attempts are in conflict with thermodynamics 
and in contradiction with experiment. This is because diffusion processes make the deviations 
from thermal equilibrium even more dramatic. 

One may conclude that the Lifshitz theory, if used within its application region, is free of 
problems and paradoxes. Problems and paradoxes arise when the basic applicability condition of 
the Lifshitz theory, thermal equilibrium, is not taken into account by using the Drude model or 
introducing the diffusion current. What is more, phenomenologically the Lifshitz theory leads 
to both thermodynamically and experimentally consistent results if the processes of electric 
conductivity connected with both drift and diffusion currents are simply disregarded. Future 
progress should show if these processes have any relation to the phenomenon of dispersion forces. 
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